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I. INTRODUCTION 



A. BACKGROUND 

Lord Rayleigh [Ref. I] was seemingly the first to investigate the mechanics of 
flows with curved streamlines. His analysis showed that inviscid curved streamline 
flows are stable providing that circulation increases with radial distance (for flow 
between concentric cylinders this condition is met when the outer cylinder spins more 
quickly than the inner cylinder). 

Later. Taylor [Ref 2] studied viscous flow between concentric rotating cylinders 
and showed that when the inner cylinder is rotated while holding the outer cylinder 
stationary, the onset of instability is predicted by the dimensionless parameter: 



where U is the tangential velocity of the inner cylinder, d is the gap between cylinders, 
n is the radius of curvature of the inner cylinder, and v is kinematic viscosity. The 
instability is manifest in the formation of regularly spaced, counter-rotating vortices 
whose axes of rotation lie in the circumferential (main flow) direction (see Figure 1.1). 
Furthermore, Taylor concluded that these vortices form in the laminar flow regime but 
do not induce turbulent main flow [Ref 3]. 

The flow of a viscous fluid in a curved channel was studied by Dean [Ref 4]. He 
developed an analytical solution for flow in a curved channel approximated by 
concentric cylinders where the radius of the inner cylinder was much greater than the 
gap between cylinders. Instead of rotating the inner cylinder Dean imposed a pressure 
gradient along the channel wall and predicted that flow instability would occur when a 
critical value of a dimensionless parameter similiar to the Taylor number was exceeded: 



Ud /d 

Ta (Taylor number) = — J — ^41.3, 



(eqn l.l) 



V 



r. 



Ud /d 

De (Dean number) = — v — ^36. 

V r. 



(eqn 1.2) 



where d represents the channel half-width. 
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Figure 1.1 Taylor Vortices Between Rotating Cylinders. 
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Gortler [Ref. 5] went on to investigate viscous flow in a curved channel and 
concluded that the flow instability predicted by Dean would take a form analagous to 
the Taylor vortices of concentric cylinder flow (see Figure 1.2). Regularly spaced 
counter rotating vortices form along axes of rotation parallel to the main flow 
direction. Yih and Sangster [Ref 6] concluded that the centrifugal force coupled with 
the velocity gradient near the concave wall caused instability due to a stratification in 
specific weight that led to vortex formation. 

The secondarx’ flows predicted by Taylor and Gortler (Taylor-Gortler vortices) 
have been confirmed experimentally using various flow visualization techniques, and 
hot wire anemometrx' [Ref 7]. The presence of these vortices in the laminar flow 
regime led to speculation of the effect of this secondarx' flow on heat transfer over 
curved surfaces. Extensive research has been done xvhich yields good experimental 
agreement with analytical solutions for laminar flow in curved channels of rectangular 
cross section [Refs. 8,9.10], and for turbulent flow in a curved rectangular channel of 
low aspect ratio (ratio of the spanwise dimension to the channel height) 
[Refs. 11,12,13]. Research in these areas indicates that there is enhancement of heat 
transfer in these flow regimes xvhen Taylor-Gortler vortices form near the concave xvall. 
However, information relating to the fluid mechanics and heat transfer in the transition 
and turbulent flow regimes for channels of high aspect ratio is verx’ limited. 



B. APPLICATIONS 

Studying heat transfer over a curved surface is just one aspect of the complex 
fluid dynamics and heat transfer problems that exist in an operating gas turbine. The 
high temperatures characteristic of a combustion engine constrain gas turbine 
performance because of metalurgical limitations of turbine blades which in turn 
restricts the elTieiency of the turbine. Information gathered in this type of 
experimentation could lead to a better understanding of the problem and facilitate the 
design of gas turbines of higher efliciency and power to weight ratio. 



C. INTENT OF THIS STUDY 

The purpose of this inx’estigation was to study the effect of curvature on heat 
transfer in a curved rectangular channel of large aspect ratio by gathering extensive 
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data relating Nusselt number to Reynolds number for flow over walls of constant heat 
flux of straight, convex and concave geometrx’. This information will provide a firm 
data base for validating computer models of the complex How' and heat transfer in 
curv’ed rectangular channels similar in geometry. 

Data for hydraulically fully developed flow in a straight rectangular channel 
heated either from above or below was taken to provide a basis for comparison with 
the flow over a heated convex or concave wall. Empirical correlations similiar to the 
familiar Dittus-Boelter correlation for Nusselt number as a function of Reynolds 
number were calculated from the data so that the heat transfer enhancement could be 
described quantitatively and compared to previous data. 

A secondaix' goal was to determine the relative eflect of buoyancy on heat 
transfer in a curved channel. The effect of buoyancy on heat transfer in forced flow 
must be considered if the effects of free and forced convection are comparable because 
the combined heat transfer is then dependent on the direction of free convection in 
comparison to that of forced convection. 
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II. EXPERIMENTAL APPARATUS AND PROCEDURE 



A detailed description of the materials and construction methods used to build 
the experimental channel and associated instrumentation is contained in G. Galyo's 
thesis [Ref. 14; pp. 25-43]. Some minor changes had to be made to the existing 
hardware and computer software and are described herein. 



A. DESCRIPTION OF THE APPARATUS 

The channel was built of plexiglass, having a channel width of 25.4 cm and a 
height of 0.635 cm, yielding an aspect ratio of 40:1. The channel was rigidly supported 
on a wooden frame and insulated over its entire length with three layers of half inch 
foam insulation. The curved section subtends ISO® of arc, and the convex wall (inner 
curved wall) had a radius of curvature of 29.7 cm (see Figure 2.1). Air was pulled 
through the channel by a compressor joined to the duct by two inch PVC pipe through 
a throttle valve. 

Four 30.48 cm x 25.4 cm wire wound rubber heaters laminated to copper plates 
were mounted on the inside walls of the channel to accomplish heating of the air. Two 
heaters were mounted in the straight test section, one on the upper wall and one on the 
lower, at a sufficient distance from the channel entrance to insure that flow was 
hydrodynamically fully developed prior to flowing over the straight heated sections. 
The two remaining heaters were installed in the curved test section, one on the concave 
wall and one on the convex wall, with the leading edges of the heaters beginning at an 
angle of 101.5° from the curved section entrance; this allowed vortices to form prior to 
flowing over the curved heated sections. The heaters were constant power density type 
capable of 600 watts each, and could be energized individually or in any combination 
of pairs by use of two power supplies. 

Thirty copper-constantan thermocouples were imbedded in each of the copper 
plates at the copper plate-rubber heater interface so that local temperatures could be 
measured (see Figure 2.2). Four thermocouples were mounted laterally across the 
channel entrance to measure inlet temperature (see Figure 2.1). Similiarly. four 
thermocouples had been mounted at the exit of the channel; however, experimentation 
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Figure 2.1 E.xperimcntal Rectangular Channel. 



IS 




Outlet 



Figure 2.2 Thermocouple Placement in the Heated Plates. 
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showed that there was insuflicient mixing of the air up to the exit which resulted in 
erratic outlet temperature measurements [Ref 14; p. 67]. Instead, eighth inch 
honeycombing w'as installed at the exit of the channel to promote mixing, and the four 
exit thermocouples were placed in the PVC pipe near the end of the channel outlet 
transition bell (see Figure 2.1). Two layers of half inch foam insulation were applied to 
the PVC to a distance of approximately fourty-eight centimeters past these 
thermocouples to minimize heat loss to the surroundings prior to measuring outlet 
temperature. One thermocouple was placed downstream of the flow measuring orifice 
and used to calculate the density of the air flowing through the orifice. Three groups 
of five thermocouples each were connected in parallel between the layers of insulation 
surrounding the heated sections so that losses to the surroundings through the 
insulation could be estimated. 

.Any combination of eighty voltages from the thermocouples and power supplies 
could be read by the Hewlett Packard 3054A automatic data acquisition system used 
for this research. Coupled to the data acquisition system were a Hewlett Packard 9S26 
computer and a printer which provided control of the data acquisition system and 
facilitated recording, reducing, and printing of data. 

A flow measuring orifice was installed in accordance with the AS.VIE Power Test 
Code [Ref 15], and water and mercur\' manometers were used to measure the pressure 
upstream of the orifice and differential pressure across the orifice so that mass flowrate 
could be calculated. 



B. EXPERIMENTAL PROCEDURE 

The procedures followed for testing were the same irrespective of the specific test 
configuration. The procedures are: 

1. Zero the manometers, shut the system throttle valve, then start the compressor. 

2. Use the throttle valve to adjust the differential pressure across the flow 
measuring orifice to that corresponding to the desired flowrate or Reynolds 
number. Determining the approximate differential pressure corresponding to a 
specific Reynolds number takes some experience and the use of the 
approximation: 
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(eqn 2.1) 



(AP)2 ^ I — V (^P), 

3. Energize the heater and adjust the voltage for the approximate value of 

- Tj^) desired. In these experiments a value of six to twelve degrees 
centigrade was used for the temperature diflerential, the objective being to 
maximize "while limiting the highest local plate temperature to a 

maximum of about 100°C. Maximizing this temperature differential lessens the 
uncertainty of the calculated Nusselt number. 

4. Experience showed that the channel reached steady state about four hours after 
first energizing the heater, and approximately one hour to steady out between 
consecutive data runs at different flowrates. When all of the following criteria 
were met, steady state was considered to have been achieved; 

a. Change in the quantity (T -T. ) of less than two percent over a ten 

minute period. ‘ 

b. Change in the heated wall temperature of less than one percent over a ten 
minute period. 

c. Change in heated wall insulation temperatures of less than one percent over 
a ten rninute period. 

5. When the above criteria for steady state are met, run the data gathering and 
reducing program (see Appendix A) several times in succession to verify that 
steady state has been reached. Three to six consecutive data collection runs at 
steady state meeting the criteria of (4) above are required to determine an 
average Nusselt number versus Reynolds number relationship. Generally the 
run-to-run variation in Nusselt number and Reynolds number were much less at 
lower flowrates. 

6. When a sufficient number of consistent data runs have been recorded, throttle 
to adjust flowrate to the next desired Reynolds number and adjust the heater 
voltage to achieve the desired plate temperature. 

The atmospheric pressure, orifice differential pressure and pressure upstream of 
the orifice had to be read from a barometer and two manometers and entered into the 
data reduction program interactively; these parameters were not read by the automatic 
data acquisition system. 
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During this study it was noted that small (one degree) ambient temperature 
perturbations had a noticeable elTect on temperatures measured throughout the test 
system; therefore, close control of ambient temperature oscillations was required so 
that temperature changes in the test system could be minimized and consistent data 
measured from run to run. 

The initial configuration of the test channel was such that buoyancy induced flow 
opposed forced flow in the curved test section. After data was taken over the full 
range of flow that the compressor was capable of, the rig was inverted and data was 
taken with buoyancy induced flow aiding forced flow. The purpose of this aspect of 
the study was to determine whether natural convection contributed significantly to the 
total heat transfer in the curved test sections. 



III. DATA REDUCTION 



Each heater plate provided a constant surface heat flux. The heater plates were 
energized by regulated D.C. power supplies through stable precision resistors of low 
resistance value (see Figure 3.1). The voltage drop across both the heater element and 
the precision resistor were read and recorded by the data acquisition system. Then the 
total power supplied by the heater was determined from the relation: 

V. X V 

Qp = ^ ^ (eqn3.1) 

^pr 

where R was a known value of resistance included as a constant in the data reduction 

pr 

program. The insulated, unheated walls were treated as adiabatic because the heat loss 
through them was negligible compared to the total heat input from the heater plate. 
Therefore, the bulk temperature varied linearly from the leading to the trailing edge of 
the heated plate. 

Thermocouple voltages were converted to centigrade temperature using individual 
thermocouple calibration curves. All temperatures were recorded then combined to 
calculate specific average temperatures desired. 

The actual heat loss to the surroundings was difficult to account for. Typically 
the estimated total heat loss by conduction through the insulation and by radiation 
was only about ten percent of the actual total heat loss. The major difficulty in 
making an accurate prediction of the heat loss was due to not being able to determine 
the area of the channel over which heat was lost to the surroundings. For the purpose 
of estimating the heat loss, the area of the heated plate was used as the effective area 
of loss; however, the plexiglass section of the channel adjacent to the heater plates acts 
as extended heat transfer surfaces, and the other copper heater surfaces in the channel, 
which have a low resistance to thermal conduction, act like heat sinks. Therefore, 
instead of assuming that the energy transferred to the air was equal to the energy input 
from the plate minus the estimated energy lost to the surroundings over the area of the 
plate, the energy actually transferred to the working fluid was calculated as follows; 
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Figure 3.1 Power Measuring Circuit. 
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(eqn 3.2) 



Q . = me AT 

^air p 

where AT = T - T.„ 

out in 

The corresponding plate average heat transfer coefficient and Nusselt number based on 
the hydraulic diameter of the duct were then calculated as follows: 



h = 



Q • 



Apl (T.h ■ T.) 



— hO.. 

air 



(eqn 3.3) 



(eqn 3.4) 



where the hydraulic diameter of the channel, is a geometric constant, 



D 



hd 



4 X Cross Sectional Area of the Duct 
Wetted Perimeter of the Duct 



(eqn 3.5) 



Data on atmospheric pressure, orifice dilTerential pressure, and orifice upstream 
pressure was used with the flow measuring orifice characteristics and air density based 
on air temperature at the orifice to calculate mass flowrate, m. Reynolds number 
based on hydraulic diameter w'as then found by the relation: 



Re 



hd 



”Ph.i 

Acl>air 



(eqn 3.6) 



where was the viscosity of the working fluid at the average bulk temperature 
measured in the test channel. Similiarly. Reynolds numbers based on the pipe diameter 
and channel height were also determined. For the concave and convex sections only, 
the dimensionless parameter. Dean number, was calculated from: 



De 



Act'd, fi 



(eqn 3.7) 
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IV. EXPERIMENTAL RESULTS 



Several sets of data were taken over the approximate range of Reynolds numbers 
from 6000 to 20,000. A computer program (see Appendix A) was used to reduce the 
raw data and calculate the quantities of interest, particularly plate average Nusselt 
numbers for corresponding Reynolds numbers, both based on the hydraulic diameter of 
the channel. Three to six separate data sets at each Reynolds number were taken, then 
averaged to give the value of the quantities of interest. 

After extensive data had been taken for the channel configured such that free 
convection opposed the direction of forced convection, the channel was inverted and 
several sets of data were taken for the convex and concave plates with the directions of 
free and forced convection coincident so that the elTect of buoyant flow could be 
determined. 

The significant parameters and corresponding uncertainty for the dilTerent 
configurations follows in Tables 1 through 4. Graphic comparisons of plate averaged 
Nusselt number versus Reynolds number are shown in Figures 4.1 through 4.7. These 
graphs also include calculated uncertainty ranges for the data, the best fit correlations 
for each configuration, data for comparison from other sources, and some of the 
correlations that have been documented in research similiar to this experimentation. 



A. HEATING THE STRAIGHT UPPER PLATE 

.A, total of sixty-nine data pairs were taken over an approximate range of 
Reynolds numbers of 6S00 to 20,000 (see Table 1). One means used to analyze the 
data was the statistics software, MINTTAB [Ref 16]. It was observed that the data at 
Reynolds numbers below 7000 exhibited a relatively large deviation from the regression 
in comparison to the data for Reynolds numbers above 7000. This result agrees with 
previous research done by J. Wilson [Ref 17; p. 50] which revealed that transition flow 
in the channel persisted up to a Reynolds number of approximately 7100; therefore, 
data at Reynolds numbers below 7000 was not included in the regression calculated for 
the upper straight plate. 
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Using a least squares regression, a correlation Tor the heat transler from the 
upper plate was determined to be: 

= 0.024 Re^/^5 (eqn 4.1) 

where Reynolds and Nusselt numbers are based on the hydraulic diameter. The 
statistical correlation of the regression was 0.998. The experimental uncertainty of the 
data w'as approximately ± 5% for Nusselt number and ± 6% for Reynolds number. 

For comparison, consider the well known Dittus-Boelter correlation based on 
hydraulic diameter and a Prandtl number of 0.71: 

= 0.02Re^/-80 (eqn 4.2) 

At first glance this small diOerence between the correlations seems insignificant; 
however, a look at Figure 4.1 shows that the Dittus-Boelter correlation overestimates 
the actual Nusselt number by about thirty to thirty-seven percent over its range of 
applicability (Re^ 10,000). 

Also shown for comparison on Figure 4.1 are interpolated values from the Kays 
and Leung [Ref 18] analytical solution for infinite parallel plates. These values too lie 
outside the uncertainty limits of the actual data, but only overestimate Nusselt number 
by ten to fifteen percent when compared to the finite flat plates of this research. 

Data previously taken by G. Galyo [Ref 14: p. 58] is included in this plot. His 
data agrees with that of this research to within the bounds of uncertainty. It must be 
noted that the data taken by Galyo at Reynolds numbers above about 17,000 is less 
certain than the other data due to the dilTiculty he encountered in recording an 
accurate channel outlet temperature at high flowrates [Ref 14: p. 67], 

B. HEATING THE STRAIGHT LOWER PLATE 

A total of sixty-four data pairs were recorded over an approximate range of 
Reynolds numbers of 6000 to 19000 (see Table 2). As in the previous case for the 
upper plate, the data recorded for Reynolds numbers less than about 7000 were not 
included in the least squares regression because of their relatively high deviation from 
the regression. The best fit regression correlation was: 
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NUj^d O.OlSRCj^d 



0.78 



(eqn -1.3) 



The statistical correlation of the regression was 0.998, and the experimental uncertainty 
of the data was essentially the same as that for the upper flat plate. 

Figure 4.2 shows that over the full range of testing there is essentially no 
difference in the results for the two flat plate cases. Data previously taken by G. 
Galyo [Ref 14: p. 59] is shown as well for comparison, and is well within the 
experimental uncertainty of the data recorded in this study except at the highest 
Reynolds number where, as mentioned previously, G. Galyo's data was highly 
uncertain. 



C. HEATING THE OUTER CURVED (CONCAVE) PLATE 

Seventy-two data pairs were taken over an approximate range of Reynolds 
numbers of 6000 to 20,000 (see Table 3). A least squares regression of the data 
resulted in the relation: 

" 0.05Rej^d°’^^ 

The statistical correlation of the regression was 0.999, and the experimental uncertainty 
was the same as that for the two previous cases. 

Figure 4.3 revealed that the previous data reported by G. Galyo [Ref 14: p. 60] 
coincides well with the data recorded in this study except, again, at the highest 
Reynolds number where G. Galyo had difficulty determining duct outlet temperature 
accurately. Also shown for comparison are extrapolated correlations from Kreith 
[Ref 12: p. 1253], and Brinich and Graham [Ref 11; p. 35]. 

In contrast to the results for the flat plates there was no significant deviation 
from the regression even at the lowest Reynolds number tested for the concave plate. 
This supports previous research by J. Wilson [Ref 17: p. 50] which indicated that the 
limiting Reynolds number for transition flow in curved sections was less than that for 
flat plates and approximately equal to 5300. 
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The data for the flat plates was used as a baseline for comparison with the 
concave plate, and heat transfer enhancement for the concave section was reflected by 
an increase in Xusselt number of about twenty-one percent at a Reynolds number of 
10,000 to fourteen percent at a Reynolds number of 20,000 (see Figure 4.4). 

There is a tendency to extrapolate this data to where the correlations for the 
concave plate and the flat plates converge in order to estimate a Reynolds number 
where the effect of curvature is no longer significant. For these correlations 
convergence occurs at a Reynolds number of about 205,000 and a Xusselt number of 
approximately 232; however, research by Brinich and Graham indicated that the 
enhancement in heat transfer in a concave section was significant even at Reynolds 
numbers as high as 250,000. 

D. HEATING THE INNER CURVED (CONVEX) PLATE 

A total of seventy-three data pairs were recorded over a range of Reynolds 
numbers of 6000 to 20,000 (see Table 4). A least squares fit of the data gave the 
relation: 

= 0.017Rej^/-8 (eqn4.5) 

at a statistical correlation of 0.997, and an e.xperimental uncertainty the same as for the 
previous cases. 

Figure 4.5 revealed that, as was the case for the concave plate, there was no 
significant deviation from the regression even at the lowest Reynolds number tested. A 
comparison of this data with that reported by G. Galyo [Ref 14: p. 61] showed that 
there is agreement to within the bounds of experimental uncertainty except at the 
highest Reynolds number. For comparison, data reported by Kreith [Ref 12: p. 1253], 
and Brinich and Graham [Ref II: p. 36] for a convex surface were included in the 
graph. 

Using the flat plate data as a reference of comparison, the reduction in heat 
transfer on the convex plate was manifest in a Xusselt number consistently six percent 
lower than that of the flat plates over the entire range of flow rates. This did not 
represent a large difference between the convex and flat plate cases; however, six 
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percent was at least one percent more than the conservatively calculated experimental 
uncertainty of the data. Therefore, this data implied that a decrease in heat transfer on 
a convex surface should be expected, though probably small. 

A comparison of the convex and concave surfaces (see Figure 4.6) pointed out 
that the concave surface heat transfer is approximately thirty percent greater than that 
for the convex surface at a Reynolds number of 8000, and about twenty percent greater 
at a Reynolds number of 20,000. 



E. EFFECT OF BUOYANCY 

Tables 3 and 4 include data taken when the test channel was oriented such that 
free convection aided forced convection. Graphical comparison for this case is shown 
in Figure 4.7. 

Data taken for this configuration was within the bounds of experimental 
uncertainty, and showed no pattern in differing from the case where the direction of 
free convection opposed that of forced convection. These results indicated that there 
was no significant buoyant effect over the full range of How rates studied. 
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6841 

7856 

8588 

9327 

10168 

10936 

11947 

12082 

12981 

13969 

14997 

15950 

17169 

17839 

18927 

20008 



TABLE 1 

SUMMARY OF STRAIGHT UPPER WALL RESULTS 









17.30 


±0.0560 


±0.0481 


19.92 

21.46 

23.06 

24.75 


±0.0576 


±0.0467 


26.18 

28.11 

28.38 

29.25 

30.96 


±0.0571 


±0.0485 


32.66 

34.62 

36.85 

37.30 


±0.0570 


±0.0495 


39.01 

40.27 


±0.0570 


±0.0509 
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TABLE 2 

SUMMARY OF STRAIGHT LOWER WALL RESULTS 







R'hd 




6071 


15.03 


±0.0581 


±0.0452 


7065 


17.85 






S116 


19.72 


±0.0573 


±0.0458 


S96S 


21.26 






9997 


23.63 






10934 


25.35 


±0.0571 


±0.0464 


12014 


27.17 






129S4 


29.11 






14047 


30.89 


±0.0571 


±0.0473 


14989 


32.67 






15981 


34.10 






16984 


35.45 


±0.0570 


±0.0479 


18035 


36.21 






19002 


37.46 






19458 


39.83 


±0.0570 


±0.0474 



TABLE 3 





SUMMARY 


OF CONCAVE WALL RESULTS 








6Rej^d'R^hd 




5884 


20.04 


±0.0582 


±0.0448 


6025* 


19.98* 






7050 


22.89 






8016 


24.73 


±0.0573 


±0.0450 


8044* 


24.80* 






9094 


26.89 






10025* 


28.08* 






10045 


28.71 






11023 


30.72 


±0.0571 


±0.0451 


11977 


32.72 






13028 


34.90 






13993 


36.19 


±0.0571 


±0.0459 


14982* 


38.05* 






15021 


37.95 






15996 


40.11 






17053 


42.51 


±0.0570 


±0.0464 


17943 


42.90 






19169 


45.26 






20189 


46.80 


±0.0571 


±0.0478 



* Indicates measurements that were taken with the test channel inverted 
such that the directions of free and forced convection were coincident. 
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TABLE 4 

SUMMARY OF CONVEX WALL RESULTS 











6065 


15.04 


±0.0581 


±0.0492 


6088-^ 


15.15* 






7050 


16.85 






8009 


18.87 


±0.0574 


±0.0498 


8034"= 


18.88* 






9058 


19.90 






10020 


22.12 






10040* 


22.48* 






11054 


24.32 


±0.0571 


±0.0525 


12003 


24.92 






12988 


27.83 






13944 


29.42 


±0.0571 


±0.0545 


14871* 


30.57* 






15003 


30.72 






15981 


32.20 






17009 


32.84 


±0.0570 


±0.0517 


17898 


35.43 






19015 


36.25 






20210 


37.73 


±0.0571 


±0.0524 



* Indicates measurements that were taken with the test channel inverted 
such that the directions of free and forced convection were coincident. 
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Figure 4.1 Straight Upper Plate Results. 
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Figure 4.2 Straight Lower Plate Results. 
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Figure 4.3 Concave Plate Results. 
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Figure 4.4 Comparison of Straight Upper and Concave Wall Results. 
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Figure 4.5 Convex Plate Results. 
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Figure 4.6 Comparison of Concave, Conve.x and Straight Plate Data. 
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Figure 4.7 Conve.K Concave Buoyancy Comparison. 
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V. DISCUSSION 



A. FLAT PLATES 

The experimental results for the upper and lower straight sections of the channel 
were essentially the same, and are not considered separately in this discussion. 

The Dittus-Boelter correlation for Nusselt number based on the hydraulic 
diameter of the test channel and a Prandtl number of 0.71 yielded Nusselt numbers 
thirty to thirty-seven percent higher than those recorded in this experimentation. The 
Dittus-Boelter correlation is an empirical appro.ximation that was developed to 
estimate local Nusselt numbers for circular tubes. The correlation itself is only an 
approximation, and when modified for non-circular ducts where convection coerTicients 
var}’ about the peripher\’, the approximation becomes less certain. Here the correlation 
was used to estimate the average Nusselt number over the length of the heated section. 
To a first approximation the Dittus-Boelter correlation can be used to estimate the 
average Nusselt number because the ratio of the plate length to the hydraulic radius 
was greater than ten; however, this assumption adds to the uncertainty of the resulting 
convection coelTicient. 

The solution for parrallel plates of infinite extent developed by Kays and Leung 
gave Nusselt numbers closer to the experimental results, but was still high by ten to 
fifteen percent, though ten to fifteen percent accuracy when estimating heat transfer 
coefficients in the turbulent flow regime is not unusual. Even though the straight 
section of the test channel approximated fiat parallel plates of infinite extent, some 
sidewall influence on flow was expected and may account for a part of the disparity of 
the comparison. The graphical comparison shown in Figure 4.1 was developed by 
linear interpolation of the data from [Ref 18: p. 552], though there was no evidence to 
support the assumption that the results of the Kays and Leung solution are a linear 
function of either Prandtl or Reynolds number; therefore, additional inaccuracy may 
have been introduced into the comparison. 
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B. CURVED PLATES 



Experimental and analytical investigation of flow and heat transfer on concave 
and convex surfaces in a rectangular duct of high aspect ratio has thus far received 
little attention; consequently, there is scant previous work available for comparison. 
Some results are available for channels of modest aspect ratios, but these do not 
generally compare well due to the significant sidewall effect introduced at low aspect 
ratios. 

Studies of the flow and heat transfer characteristics in high aspect ratio channels 
has been done at the Naval Postgraduate School [Refs. 9,14,17,19,20,21,22]. Results of 
this previous research were comparable to the results of this experimentation and were 
highly depended upon to validate these results. Only Galyo's work is shown for 
comparison in this report, but his results agreed within the limits of experimental 
uncertainty with those of Wilson [Ref 17], and Daughety [Ref 22]. The major 
improvements of this experimentation over that previously done were the decreased 
experimental uncertainty, especially at the highest Reynolds number, and the 
preponderance of data that permitted empirical correlations to be determined to a very 
high level of confidence. 

In 1955 Kreith [Ref 12: p. 1253] investigated the influence of curvature on heat 
transfer in a rectangular channel of aspect ratio equal to three. For a Prandtl number 
of 0.71, graphical comparison with the results recorded in this experimentation showed 
that Kreith's results were approximately thirty percent lower (see Figures 4.3 and 4.5); 
however, Kreith's test apparatus not only had a much lower aspect ratio, it was also 
heated by a heating element that was centered between the upper and lower walls of 
the channel instead of heating the walls themselves. 

Research done in 1977 with a channel of aspect ratio six by Brinich and Graham 
[Ref II: pp. 35-36] was shown for comparison in Figures 4.3 and 4.5. The range of 
Reynolds numbers over which they performed their experiments was approximately 
S4,000 to 253.000. For the concave surface. Brinich's and Graham's data extrapolated 
to a Reynolds number of 15,000 yielded a N'usselt number about fifty percent higher 
than that recorded in this experiment. Likewise, an extrapolation to a Reynolds 
number of 15.000 for the convex plate showed that there was about a thirty percent 
difference in the results of the two experiments. These differences are rather large, but. 
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again, the two studies differed greatly in that the aspect ratios and the range of 
Reynolds numbers investigated were much different. 

Similiar research was done in 1978 by Uchida and Koizumi [Ref 13: pp. 
1714-1716], with channels of aspect ratio up to four; however, the full report has not 
been translated into English. From figures, graphs, and equations included in their 
report it appears that their test channel was heated along its entire periphery’ and that 
the published results were in the form of an average Nusselt number for the 
combination of a convex and concave wall. 

It is unlikely that an average heat transfer coefficient or Nusselt number can be 
determined when both the convex and concave walls are heated by simple 
superposition of the individual correlations for the convex and concave walls because 
of the complex nature of the flow in a curved channel. Therefore, there was no 
manner in which the results detailed in Ref 13 could justifiably be compared to the 
results of this experimentation. 

Even though the test equipment was such that a pair of heaters could be 
energized to supply a constant heat flux on opposite walls of the channel, no 
conclusive experiments could be performed to determine an average Nusselt number 
for this configuration. 

Referring to (eqn. 3.3) and denoting hj as the average heat transfer coefficient for 
single plate heating, and h^ as the average heat transfer coefficient when heating a pair 
of opposite walls, then: 




K\-o ^ 2 



(eqn 5.1) 



where represents the average temperature when one wall (inner or outer) is heated, 
and the heat flux from both the inner and outer plates was assumed equal when 
heating a pair of plates. The average wall temperature term in the numerator of (eqn. 
5.1) can be recast as: 



(T„, 



+ T„,) 2 - T„,, + 5T 



(eqn 5.2) 
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where 6T represents the difierence between the average temperature when heating a 
pair of plates and the single plate average temperature. Combining (eqn. 5.1) and (eqn. 
5.2) yields: 




+ 6T) 

(^vvh ■ 

1 + - T„,) 

1 

1 + 6T/AT 

i 



(eqn 5.3) 



(Eqn. 5.3) reveals that only when the temperatures of both plates are the same, such 
that St = 0, when heating a pair of plates, can a true average heat transfer coefTicient 
for both walls be obtained. Strictly speaking the temperature at all corresponding 
points opposite each other on the plates should be equal to yield an average heat 
transfer coefTicient for the combination of both walls. The heaters used in this 
experimentation were constant heat flux heaters, and temperatures at corresponding 
points on opposite walls could not be matched; therefore, conclusive data for the case 
of heating a pair of plates could not be gathered using the existing equipment. 

Though local heat transfer coefficients and Nusselt numbers were calculated at 
the leading edge, at mid-plate, and at the trailing edge of the heated plates, no 
conclusive pattern was evident in the data. The same was true for local temperatures. 
More conclusive results could not be obtained because there were too few 
thermocouples spaced too far apart; that is. the spanwise thermocouple spacing was 
greater than the diameter of the Taylor-Gortler vortices that formed. 
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VI. CONCLUSION 



The results of this experimentation, detailed in Chapter IV, indicated that 
secondary' flow in a curved rectangular channel in the form of Taylor-Gortler vortices 
influenced the heat transfer on the concave and convex surfaces over a range of flow 
corresponding to Reynolds numbers from about 6000 to 20000. 

Correlations of Xusselt number as a function of Reynolds number for the two 
flat plates, the convex plate, and the concave plate were calculated to a high level of 
certainty and showed that heat transfer was enhanced by approximately twenty percent 
on the concave surface as compared to the flat surfaces. Similiarly, heat transfer over 
the convex surface was reduced by approximately six percent in comparison to the flat 
surfaces. Furthermore, the enhancement or reduction in heat transfer on the curved 
plates was noted to be greatest at the lowest Reynolds number tested, and tapered ofT 
with increased flow rate. 

Comparisons made of the results obtained in this experimentation to that of 
other pertinent research supported the conclusion that heat transfer is enhanced on a 
concave surface. Experiments performed for the purpose of determining the relative 
significance of buoyancy induced flow indicated that even at the lowest Reynolds 
number of 6000, buoyant free convection was negligible compared to forced 
convection. 
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VII. RECOMMENDATIONS 



A. EQUIPMENT 

1. Small oscillations in ambient temperature ( ± T) noticeably alTected 
temperatures in the channel and contributed to the difierences in the results in 
successive data taking runs. Maintaining close control over ambient 
temperature would make experimental results more reproducible and decrease 
the time necessary to gather data. 

2. Existing equipment allowed for testing only up to a Reynolds number of about 
20000, though the corresponding Reynolds numbers for flow through a gas 
turbine are much higher. A larger test channel and higher capacity air 
compressor would be required for experimentation at higher flow rates. 

3. A great amount of time was spent waiting approximately four hours for the test 
channel to reach steady state after first energizing the heater. Installing guard 
heaters would significantly reduce the warm up time to steady state. 
Additionally, guard heaters would prevent heat loss through the insulation 
surrounding the heaters. 

4. Pressure sensors that could be read directly by the automatic data acquisition 
system several times during the course of a data run would increase the 
accuracy of flow rate measurements and permit successive data runs to be taken 
without the presence of an operator. 

5. More thermocouples should be installed closer together to improve the accuracy 
of temperature measurements and to facilitate studying the pattern of local 
variations in temperature and heat transfer in the spanwise and flow-wise 
directions. 

6. The single most significant contributor to experimental uncertainty in the 
calculated average Nusselt number and hydraulic Reynolds number was the 
uncertainty in the channel cross section. The uncertainty in experimental 
results could be greatly reduced if a channel were constructed paying special 
attention to the cross section dimensions. 
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B. FURTHER RESEARCH 

1. Experimentation similiar to this but at higher flow rates should be done so that 
a large data base could be established and heat transfer correlations could be 
calculated for a wide range of application. 

2. Experiments at lower flow rates where buoyancy is significant should be done 
so that a suitable length parameter can be deduced for calculating the Grashof 
number for curved channel flow. 

3. Different curvatures should be studied to determine the effect of the degree of 
curvature on heat transfer in a curved channel. 

4. Thermostatic control should be included in the experimental equipment that 
could maintain equal plate average temperatures when heating a pair of heaters 
on opposite walls so that experimental results for channel average heat transfer 
coeflicients could be gathered. 
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APPENDIX A 

DATA REDUCTION PROGRAM 



! FILE NOME: MAIN 
! LAST REVISION DATE: 19 OCT 1986 

! EDITED BY: LT. J.R. HAWK III 

! PROGRAM FOR GATHERING AND REDUCING DATA WHEN HEATING ONE COPPER PLATE, 

! WITH B MINUTES BETWEEN RUNS. PROGRAM IS MODIFIED TO ACCOUNT FOR 
! BROKEN THERMOCOUPLE CHANNELS B, 16, 21, AND 26 WHICH EFFECT OUTER 
I CURVED HEATER, THERMOCOUPLE CHANNEL 07 WHICH EFFECTS INNER CURVED 
! HEATER, AND THERMOCOUPLE CHANNEL 11 WHICH EFFECTS OUTER STRAIGHT HEATER 

COM /Co/ D(7) ,Aa(76,2) 

DIM Emf (79) ,T( 79) 

! CORRELATION FACTORS FOR EMF TO DEGREES CELSIUS FOR CU-CO T/C. 

DATA 0. 10086091 ,25727 .90369, -76730B.829B,7802B59B. 81 
DATA -9207086589,6. 97683E>1 1 , -2 . 66 1 92E > 1 3 , 3 . 900 78E H 0 
READ D(0 
PRINTER IS 701 
BEEP 

INPUT "ENTER MONTH, DAY, AND TIME ( MM : DD : HH : MM : SS ) * , T ime$ 

OUTPUT 709;"TD";Time$ 

BEEP 

INPUT "ENTER INPUT MODE ( 0 = 3050A , 1 »F ILE ) Mm 

IF Im»l THEN 

BEEP 

INPUT "ENTER NAME OF EXISTING DATA FILE " ,01df i le$ 

PRINT USING "20X, ""These Results Are From Data File: " " , 1 OA" ; Oldf i le$ 

PRINT 

PRINT 

ELSE 

BEEP 

INPUT "ENTER NAME OF NEW DATA FILE'>" ,Neuf i le$ 

PRINT USING "20X,""Data For This Run Is Stored In Data File: "".10A";Neuf 

PRINT 
PRINT 
END IF 

! ENTER THE CHANNEL CONFIGURATION. 

PRINT USING "9X,""THE FOLLOWING CHANNEL CONFIGURATION WAS SELECTED:""" 

PRINT 

BEEP 

INPUT "SELECT HEATER TYPE (0 = STRAIGHT, 1 =CURVED ) " , I type 

IE ltype=0 THEN PRINT USING " 1 2X , " " Mea t i ng STRAIGHT Test Section.""" 

IE Itype=l THEN PRINT USING " 1 2X , " "Hea t i ng CURVED Test Section. 

BEEP 

INPUT "SELECT HEATER POSITION (O-OUTER, 1 = INNER )", Ipos 
IE lpos=0 THEN PRINT USING " 1 2X , " "Hea t i ng OUTER Plate.""" 

IE Ipos=1 THEN PRINT USING " 1 2X "Hea t i ng INNER Plate.""" 

PRINT 
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1S20 
1S30! 

ISSO 
15G0 
1S70 
1S80 
1S90 
IGOO! 

IGIO 
1G20! 

1630 
1G40 
1650 
1G60 
1670 
1680 
1690 
1700 
1710 
1720 
1730 
1740! 
1750 
1760! 
1770 
1780 
1790 
1800 
1810 
1820 
1830 
1840 
1850 
1860 
1870 
1880 
1890! 
1900 
1910 ! 
1920 
1930 
1940 
1950 
1960 
1970 
1980 
1990 
ad ) 
2000 
da) 

20 1 0 

s) 

2020 

s) 

2030 
2040 
2050 
2060 ! 
2070 
2080! 



EmERS FILES «.H I/C CRLIBRRHOM COEFFICIENTS BRSEO OH HESTER TYRE. 



s fsifc!; fFiir. is 

ENTER ^rluc^ftad.O) ,8a(I.n .8a(1.2) 

NEXT I 

assign ‘SFilec TO » 

! INSERT DATA STORAGE DISK. 

input "CHANGE DISK AND 

IF Im=l THEN ASSIGN L«File TO OldfileS 

IF Iro'O THEN , , . 

CREATE BOAT NeufiIeS.40 
assign sFile TO NeutileS 
END IF 

enter 709; Times 

lllHO -20X.--Ho„tH. 0„. .«OTI~: 

PRINT 

, for DSTR RON RESD SMO STORE OLL R«H EHF VRLUES. 

Sim/os-a® »f“» "",rFL.- 

OUTPUT 722;"F1 R1 T1 Z1 FLl 
for I’O to 79 
OUTPUT 709; "AS" 
enter 722;Emf(I) 

OUTPUT iJF i le;Emf (*) 

ENTER i^F ile;Emf (*' 

END IF 

™'cONVErT’sN°D CRLIBRHTE RLE I/C REHOINCS TO DECREES CELSIUS. 
PRINT USIHC -9X.--IHE FOLLOHING DSTP HERE RECORDED: 

PRINT 

?l■HiU^o^“4B5■J'R?^pi;rr..//TFHEN 



(U)"““;Emf(38'-Ia 
, Mip) - ““.20. 30.”" ;Emf(39da 

print USING "12X. ""Resistor Voltage 2.01729 (Ohm 

,E IFO.0 THEN PRINT USING ■ ~ -Rt.. 

IF Iposd THEN PRINT USING "12X." Precision 

PRINT 
END IF 
ELSE 



IF T/C OUT OF RANGE PRINT A OUE . 
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2090 Aemt-Emf(I) 

2100 IF Aemf<7.09E-^ OR Aemf >3 . 630E-3 THEN 
2110 T(I)-999.99 

2120 ELSE 

2130 CALL Tvsv(Emf(I),Tt) 

2M0 IF 1-79 THEN 
2150 T(I)-Tt 
2150 GOTO 2210 
2170 END IF 
2180 T(I)-FNTcorr(Tt,I) 

2190 END IF 
2200 END IF 
2210 NEXT I 
2220 ! 

2230 } PRINT ALL T/C READINGS. 

22 ^ 0 ! 

2250 PRINT 
2260 PRINT 
6 " 

2270 PRINT 
(5) 

2280 PRINT 
12 

2290 PRINT 

T(in 

2300 PRINT 
13 

2310 PRINT 
16) .T(17) 

2320 PRINT 
2 ^ 

2330 PRINT 
22) ,T(23) 

23^0 PRINT 
3C " * ” 

2350 PRINT 
28) ,T(29) 

2360 PRINT 
2370! 

2380 ! CALCULATE AVERAGE TEMPERATURE OF OUTER HEATED WALL. 

2390! 

2400 Sum-0, 

2410 IF Itype-0 THEN 

2420 FOR 1=0 TO 29 

2421 IF 1=11 THEN 2440 
2430 Sum=Sum+T(I) 

2440 NEXT I 

2450 Tuo=Sum/29 
2460! 

2470 ! CORRECTION FOR BROKEN T/C CHANNELS 5, 16, 21, AND 26 IN OUTER CURVED HE 

ATER. 

2480! 

2490 ELSE 

2500 FOR 1=0 TO 29 

2510 IF 1=5 OR 1=16 OR 1=21 OR 1-26 THEN 2530 
2520 Sum = Sum + T ( I ) 

2530 NEXT I 
2540 Tuo=Sum/26 
2550 END IF 

2560 PRINT USING ** 1 2X "OUTER INSULATION:"-" 



USING "12X, “"OUTER PLATE: 
USING “14X.“"T/C Number: 


II II M 

1 


USING "HX, ""Temp (C): 


““,6 (3D 


USING "HX,""T/C Number: 


7 


USING "14X,“"Temp (C): 


,6(3D 


USING “MX, ""T/C Number: 


13 


USING "MX. ""Temp (C): 


" ",G(3D 


USING "MX. ""T/C Number: 


19 


USING “MX. ""Temp (C) : 


"-,6 (3D 


USING "MX. ““T/C Number: 


25 


USING “MX, ““Temp (C) : 


"-,6(3D 



2 3 4 5 

.DD,2X)-;T(0),T(1 ),T(2) ,T(3),T(4),T 
8 9 10 11 

.DD,2X)-;T(6),T(7),T(8),T(9).T(10), 
14 15 16 17 

.DD,2X)-;T(12),T(13),T(14),T(15),T( 
20 21 22 23 

.DD,2X) -;T( 18) ,T( 19) ,T(20) ,T(21 ) ,T( 
26 27 28 29 

.DD,2X) -;T(24) ,T(25) ,T(26) ,T(27) ,T( 
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2B70 PRINT USING Number: 31 (Tinsol) 32 iJinsoZ) 33 

2B80 PRINT USING “l^X/’-Temp (C): - - , 3 ( 3D . DD . 3X ) - ; T ( 30 ) , T ( 3 1 ) . T 

2B90 Tinso1»T(30) 

2600 Tinso2»T(31) 

2610 Tinso3-T(32) 

2620 PRINT 

2630 PRINT USING ” 1 2X "ORIFICE TEMP (Torf) - "",30.00,"" ( C) ; T ( 33 ) 
26^0 PRINT 

26B0 PRINT USING " 1 2X INLET TEMPERATURE:""" 

2660 PRINT USING "1^X,""T/C Number: 3B 36 37 38""" 

2670 PRINT USING "MX, ""Temp (C): " " , ^ ( 30 . DO . 2X ) " ; T( 34 ) , T ( 3B ) , T ( 36 ) 

2680 PRINT 
2690! 

2700 ! CALCULATE AVERAGE INLET TEMPERATURE. 

2710! 

2720 Sum*0. 

2730 FOR M34 TO 37 
2740 Sum“Sum + T ( I 
27B0 NEXT I 
2760 T i n“Sum/4 
2770 PRINT USING 
2780 PRINT USING 
46 

2790 PRINT USING 
44) .T(4B) 

2800 PRINT USING 
B2 

2810 PRINT USING 
BO).T(BI) 

2820 PRINT USING 
B3 

2830 PRINT USING 
B6) ,T(B7) 

2840 PRINT USING 
64 

28B0 PRINT USING 
62) ,T(63) 

2860 PRINT USING 

70 n n M 

2870 PRINT USING 
68) .T(69) 

2880 PRINT 
2890! 

2900 ! CALCULATE AVERAGE TEMPERATURE OF INNER HEATED WALL. 

2910! 

2920 Sum*0. 

2930 IF Itype’O THEN 
2940 FOR IMO TO 69 
29BG Sum = Sum T ( I ) 

2960 NEXT I 
2970 Tiji=Sum/30 
2980! 

2990 ! CORRECTION FOR BROKEN T/C CHANNEL 47 IN INNER CURVED HEATER. 

3000! 

3010 ELSE 

3020 FOR IMO TO 69 
3030 IF M47 THEN 30B0 
3040 Sum = 5um + T ( I ) 

30B0 NEXT I 



"12X,"" INNER PLATE: 
"14X,""T/C Number: 

"14X,""Temp (C): 

"14X,""T/C Number: 

"14X,""Temp (C): 

"14X,""T/C Number: 

"14X,""Temp (C) : 

"14X,""T/C Number: 

"14X,""Temp (C): 

"14X,""T/C Number: 

"14X,""Temp (C) : 



41 42 

"",6(3D.DD,2X)" 
47 48 

"",6(30. DO, 2X)" 
B3 B4 

"",6(30. DO, 2X)" 
B9 60 

"",6(30. DO, 2X)" 
6B 66 

"",6(30. DO, 2X)" 



43 44 

T(40) ,T(41 ) , T(42) 
49 BO 

T(46) ,T(47) ,T(48) 
BB B6 

T(B2) ,T(B3) ,T(B4) 
61 62 

T(B8) ,T(B9) ,T(60) 
67 68 

T(64) ,T(6B) ,T(66) 



( T 1 nso3 ) " 
(32) 



,T(37) 



4B 

,T(43),T( 

B1 

,T(49) ,T( 
B7 

,T(BB) ,T( 
63 

,T(61 ) ,T( 
69 

,T(67) ,T( 
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3060 Tui»Sum/29 
3070 END IF 

3080 PRINT USING “ 1 2X . INNER INSULATION:*** 

3090 PRINT USING "1AX.“"T/C Number: 71 (TinslD 72 (Tinsi2) 73 (Tinsi3) 

II m 

3100 PRINT USING ”1^X,”-Temp (C): ” 3 ( 3D . DD , 8X ) - ; T ( 70 ) , T ( 7 1 ) , T ( 72 ) 

3110 Tinsil-T(70) 

3120 Tinsi2-T(71) 

3130 Tinsi3-T(72) 

31^0 PRINT 

31S0 PRINT USING “ 1 2X "OUTLET TEMPERATURE:""" 

3160 PRINT USING "MX,""T/C Number: 7^ 75 76 77""" 

3170 PRINT USING "MX/’"Temp (C): " " , ^ ( 3D . DD , 2X ) " ; T ( 73 ) J ( 74 ) , T ( 75 ) . T ( 76 ) 

3180 PRINT USING 

3190! 

3200 I CALCULATE AVERAGE OUTLET TEMPERATURE. 

3210! 

3220 Sum-0. 

3230 FOR 1-73 TO 76 
3240 Sum-Sum+T(I) 

3250 NEXT I 
3260 Tout-Sum/4 
3270! 

3280 I CALCULATE DIFFERENCE BETWEEN CHANNEL INLET AND OUTLET TEMPERATURES. 
3290! 

3300 Tdif f-Tout-Tin 
3310! 

3320 ! CALCULATE CHANNEL BULK TEMPERATURE. 

3330 ! 

3340 TbIk-(Tin+Tout)/2 
3350! 

3360 ! PRINT CALCULATED TEMPERATURES. 

3370! 



3380 


PRINT 


USING 


"9X,""THE FOLLOWING TEMPERATURES WERE CALCULATED:""" 






3390 


PRINT 












3400 


PRINT 


USING 


"12X, ""Average 


Outer Wail Temperature (Two) 


» ** - 


,3D 


DD, " " 


( C ) " " 


" :Tuio 










3410 


PRINT 


USING 


"12X, ""Average 


Inner Wail Temperature (Twi) 


m •* - 


,3D 


DD. " " 


( C ) " " 


" ; Ty i 








3420 


PRINT 


USING 


" 1 2X , " "Average 


Outlet Temperature (Tout) 


« - 


,3D 


DD."" 


( C ) " " 


" : I o u t 










3430 


PRINT 


USING 


" 1 2X , " "Average 


Inlet Temperature (Tin) 


m - 


,3D 


DD," " 


( C ) " " 


";Tin 








3440* 


PRINT 


• USING 


"12X, ""Channel 


Inlet and Outlet Temp Diff (Tdiff) 


» *• « 


.3D 


DD." " 


( C ) " " 


" ; T d 1 f f 










3450 


PRINT 


USING 


"12X." "Average 


Bulk Temperature (Tbik) 


m 


.3D 



DD,"" (C) ;Tblk 

3460! 

3470 ! CALCULATE DELTA T (Tdei) FOR USE IN CALCULATING H LATER. 

3480! 

3490 Tdelo=T(jo-Tbik 

3500 Tde i i = Tw i -Tb He 

3510 IF lpos=0 THEN Tdel-Tdelo 

3520 IF Ipos=1 THEN Tdel=Tdeli 

3530 PRINT USING “ 1 2X . " "Mean Temperature Difference (Tdei) - 

DD."" (C) ;Tdel 

3540! 

3550 ! CALCULATE DIFFERENCE BETWEEN OUTER AND INNER WALL. 

3560! 

3570 Tud 1 f f * ABS ( TujO' Tui 1 ) 
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3S80 PRINT USING ” 1 2X , " “Ou ter and Inner Wall Temp Difference (Tuidiff) « "\3D. 

DD,“" (C) ;Tujdiff 

3590! 

3600 ! CALCULATE LOCAL CHANNEL TEMPERATURES AT INLET, MIDDLE AND OUTLET. 

3610! 

3620 Tc 1 n- T 1 n+ ( . 5* Td 1 f f / 1 2 ) 

3630 Tcmid-Tin>(^ .S^Tdiff /12) 

3640 lcout-Tin+( 1 1 ,5^Tdiff/12) 

3650 PRINT USING '* 1 2X "Average Local Channel Temp at Inlet (Tcin) - “",3D. 

DD."" (C)""";Tcin 

3660 PRINT USING " 1 2X "Average Local Channel Temp at Middle (Tcmid) ■ ■*,3D, 
DD,"“ (C)"“";Tcmid 

3670 PRINT USING " 1 2X "Average Local Channel Temp at Outlet (Tcout) • ""tSD. 

DD , ( C ) " ” ” ; Tcou t 

3680! 

3690 ! CALCULATE AVERAGE LOCAL HEATED WALL TEMPERATURE AT INLET. 

3700! 

3710 Sum-0. 

3720 Iadd-0 

3730 IF Ipos-I THEN Iadd-40 
3740! 

3750 ! CORRECTION FOR BROKEN T/C CHANNELS 5, AND 47 CURVED HEATER ONLY. 

3760! 

3770 IF Itype-1 THEN 
3780 FOR I-O TO 6 

3790 IF (Ipos-0 AND 1-5) OR (Ipos-1 AND 1-6) THEN 3810 
3800 Sum-5um+ T ( I + ladd ) 

3810 NEXT I 
3820 Tuin-Sum/6 
3830 ELSE 
3840 FOR I-O TO 7 

3850 IF (Ipos-0 AND 1-7) OR (Ipos-1 AND 1-6) THEN 3870 
3860 Sum-Sum+T ( I + ladd ) 

3870 NEXT I 
3880 Tuin-Sum/7 
3890 END IF 

3900 PRINT USING "12X,""Avg Local Heated Wall Temp at Inlet (Twin) - "",3D. 

DD,"" (C) ;Tuin 

3910! 

3920 ! CALCULATE AVERAGE LOCAL HEATED WALL TEMPERATURE NEAR MIDDLE. 

3930! 

3940 Sum-0. 

3950! 

3960 ! CORRECTION FOR BROKEN T/C CHANNEL 16 CURVED HEATER ONLY, AND T/C 

396HCHANNEL 11 STRAIGHT HEATER ONLY. 

3970! 

3980 IF Itvpe-t THEN 
3990 FOR 1-9 TO 16 

4000 IF (lpos=0 AND I-tO) OR dpos-O AND 1=16) OR (Ipos-1 AND 1=9) THEN 4020 
4010 Sum=Sum+ T ( I ^ ladd ) 

4020 NEXT I 

4030 IF Ipos-0 THEN Turn id-Sum/6 . 0 

4031 IF Ipos-1 THEN Turn id=Sum/ 7 . 0 
4040 ELSE 

4050 FOR 1=9 TO 16 

4060 IF (lpos=0 AND I-IO) OR (Ipos-1 AND 1=9) OR (Ipos-0 AND I-ll) THEN 4080 
4070 Sum = Sun + T ( I ladd ) 

4080 NEXT I 

4081 IF Ipos-0 THEN Turn id=Sum/6 . 0 
4090 IF Ipos-1 THEN TufTud = Sum/7 . 0 
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4100 END IF 

4110 PRINT USING “12X,‘*"Avg Local Heated Wall Temp Near Middle (Twmid) ■ "",3D 

DD,”*' (C) ;Tumid 

4120! 

4130 ! CALCULATE AVERAGE LOCAL HEATED WALL TEMPERATURE AT OUTLET. 

4140! 

4150 Sum-0. 

4160! 

4170 I CORRECTION FOR BROKEN T/C CHANNEL 26 CURVED HEATER ONLY. 

4180! 

4190 IF Itype-1 AND Ipos-0 THEN 
4200 FOR 1-23 TO 29 
4210 IF 1-26 THEN 4230 
4220 Sum-Sum^'T ( I ) 

4230 NEXT I 
4240 Tuout-Sum/6 
4250 ELSE 

4260 FOR 1-23 TO 29 
4270 Sum-Sum+T( I+Iadd) 

4280 NEXT I 
4290 Tuout-Sum/7 
4300 END IF 

4310 PRINT USING •12X.**"Avg Local Heated Wall Temp at Outlet (Tyout) - “"»3D 

DD."** (C)"“‘';Tuout 

4320 PRINT 

4330 PRINT 

43401 

4350 ! ENTER ORIFICE PRESSURE DATA. 

4360! 

4370 BEEP 

4380 INPUT "ENTER PATM (InHg), DPM(inH20). PlM(inH20), REl “ ,Patm,Dpm,Plm,Rel 
4390 PRINT USING "9X,"-THE FOLLOWING ORIFICE DATA WERE ENTERED:"'*" 

4400 PRINT 

4410 PRINT USING ** 1 2X "Pa tm( i nHg ) DPM(inH20) PlM(inH20) REl""" 

4420 PRINT USING " 1 4X , 2 ( 2D . DD J OX ) , ( 2D . DD . 7X ) , 5D . D" ; Pa tm . Dpm , PI m , Rel 
4430 PRINT 
4440! 

4450 ! CONVERT PRESSURE READINGS TO SI UNITS AND PIM TO ABSOLUTE. 

4460! 

4470 Pa tm-Pa tm*3376 . 8 
4480 Pdel-Dpm*248.84 

4490 Pl=Patm-(Plm>248.84) 

4500 Torf-T(33) . 

4510! 

4520 ! PHYSICAL PROPERTIES AND CONSTANTS. 

4530! 

4540 R-286.987 !FOR AIR. 

4550 Rho=Pl/(R*(Tort>273. 15) ) !DENSITY OF AIR AT ORIF PLATE. 

4560 Gamma=l .40 ! FOR AIR. 

4570 Cp=1006 !Cp GOOD FOR Tblk BETWEEN 12 AND 33 (DEG C) . 

4580 Gc=1.0 

4590 Sigma=5.669E-8 

4600 Ecu-. 12 !EMISSIVITY OF COPPER. 

4610! 

4620 ! DIMENSIONS OF CHANNEL. 

4630! 

4640 Dc-. 006350 
4650 Wid-.2540 

4660 R1-.297 

4670 Pye t-2*- ( Dc+W id ) 

4680 Ac*Dc'^Wid 



55 



^690 Dhd=^*Ac/Pue t 
A700 Apl-. 077^192 
^710! 

4720 ! PROPERTIES OF INSULATION. 

4730! 

4740 Xins-.0127 

4750 Kins*3.8E-2 

4760! 

4770 ! ENTER ORIFICE CONFIGURATION. 

4780! 

4790 BEEP 

4800 INPUT “SELECT DIAMETER OF ORIFICE (0-.5334, 1*1.0755 ( inches) )“ ,Si le 
4810 IF Size*0 THEN Dorf*. 013548 
4820 IF Size*1 THEN Dorf*. 027318 
4830! 

4840 ! DIMENSIONS OF PIPE AND ORIFICE PLATE. 

4850! 

4860 A»(PI^Dorf‘ 2)/4 
4870 Dpipe*. 051772 
4880 Ap ipe"PI*Dp ipe ‘2/4 
4890 Beta*Dor f /Dp ipe 

4900 PRINT USING “ 1 2X . “ “Do rf ( m ) A(m*2) Beta“““ 

4910 PRINT USING “ 1 2X , ( Z . 6D , 6X ) , ( Z . 3DE , 6X ) , Z . 4D“ ; Do r f , A , Be ta 
4920 PRINT 
4930 PRINT 
4940! 

4950 ! CORRELATION FOR EXPANSION FACTOR BASED ON ID AND 1/2D TAPS. 

4960! 

4970 Y-1-( .41^.35-*Beta‘4)^(Pdel/( Gamma^P 1 ) ) 

4980! 

4990 ! CALCULATION OF FLOW COEFFICIENT (K) IN SI UNITS. 

5000! 

5010 B-. 0002+2. 794E-5/DpipeM .0038+1 .0l6E-5/Dpipe )»1 Beta ‘2>( 1 6 . 5+ 1 . 968504E+2*Dp 

ipe ) »Be ta ‘16) 

5020! 

5030 ! K1 AND K2 ARE DUMMY VARIABLES TO CALCULATE Ko. 

5040! 

5050 K1-.6014-5.3974064E"3/Dpipe‘ .25 

5060 K2=( .376+2.8971 138E-2/Dpipe‘ .25) ♦( 1 .6129E-7/( Dpipe*2+Beta ‘2+6.35E-5*Dpipe) 

+Beta*4+1 .5»Beta*16) 

5070 Ko»K1+K2 

5080 K*Ko+l000»B/Re1 ‘ .5 

50901 

5100 ! Mu CORRELATION GOOD FOR Torf BETWEEN 17 AND 41 DEGREES CELSIUS. 

5110! 

5120 Mu=4.6971429E-8+Torf+l .7194722E-5 

5130! 

5140 ! CALCULATE Mdot AND REpipe AND COMPARE TO PREDICTED RE1. 

5150! 

5160 Mdot = K + A^Y'<'(2'*Gc^Rho + Pdel ) * .5 
51 70 Repipe=(Mdot+Dpipe)/(Mu*-Apipe) 

5180 Dif ABS(Re1-Repipe)/Rel ) + 100 

5190 Re 1 =Rep ipe 

5200 IF Diff>.00l THEN 5080 

5210 PRINT USING “SX/’-THE FOLLOWING DATA WERE CALCULATED “ 

5220 PRINT 

5230 PRINT USING “ 1 2X . “ “Or i f ice Expansion Factor (Y) * ““,Z.4D.“;Y 

5240 PRINT USING “ 1 2X . “ “Or i f ice Flow Coefficient (K) * ““.Z.4D/‘:K 

5250 PRINT USING “ 1 2X . “ “Dens i ty Based on Torf (Rho) * ““,Z.4D,““ (kg/m‘ 

3 ) “ “ “ ;Rho 
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5260 PRINT USING “12X 
;Mu 



"“Viscosity Based on Torf 



5270! 

5280 

5290! 

5300 

5310 

5320 

5330! 

53^0 

5350! 

5360 

5370 

5380 

5390 

5400 

5410! 

5420 



! Mu CORRELATION GOOD FOR Tblk BETWEEN 17 

Mu-4.6971 429E~8^Tblk+1 .7194722E-5 
Red - ( Mdo t*Dc ) / ( Mu^<-Ac ) 
Rehd-(Mdot^Dhd)/(Mu*Ac) 

! CALCULATION OF POWER INTO HEATER PLATE. 
Iadd-0 

IF Ipos-1 THEN Iadd-39 
IF Ipos-0 THEN Rpr-2. 01729 
IF Ipos-1 THEN Rpr-2. 00839 
Qp-(Emf (38-‘-Iadd)*-Emf (39 + Iadd) )/Rpr 

! CALCULATION OF HEAT INPUT TO AIR. 



5430! 



(Mu) 



AND 41 



- "",2Z.3DE,“" (kg/m.s) 



DEGREES CELSIUS. 



5440 Qa 1 r-Mdot^^Cp* ( Tout~T i n ) 

5450! 

5460 1 CALCULATE LOCAL HEAT TRANSFER COEF. AND NUSSELT NUMBER. 

5470! 

5480 Ka 1 r- 7. 7257 143E"5*Tcin+. 024 165836 
5490 H i n-Qa i r/( Ap 1 «■ ( Tu i n-Tci n ) ) 

5500 Nu i n- ( H i n^^Dhd ) /Ka i r 

5510 Ka i r= 7. 7257 143E-5* Tern id+. 024 165836 

5520 Hmid-Qa i r/ ( Apl *( Tumid-Tcmid) ) 

5530 Numid-(Hmid^eDhd)/Kair 

5540 Ka i r - 7 . 7257 1 43E~5^ Tcou t + . 024 1 65836 

5550 Hou t-Qa i r/( Ap 1 *■ ( Tuou t - Tcou t ) ) 

5560 Nuout- (Hout^^Dhd ) /Ka 1 r 
5570! 

5580 ! Kair CORRELATION GOOD FOR Tblk BETWEEN 17 AND 41 DEGREES CELSIUS. 

5590! 

5600 Ka 1 r- 7 . 7257 1 43E-5^ Tb Ik + . 024 1 65836 
5610! 

5620 ! CALCULATE AVERAGE HEAT TRANSFER COEF. AND NUSSELT NUMBER. 

5630! 

5640 Havg-Qa 1 r/ ( Ap 1 *Tde 1 ) 

5650 Nuavg- ( Havg*Dhd ) /Ka 1 r 
5660! 

5670 ! CALCULATE HEAT LOSSES. 

5680! 

5690 Qlol-(Kins^Apl«-(Tinsol-Ti nso2 ) ) /X i ns 
5700 Qlo2-(Kins*Apl*(Tinso2-Ti nso3 ) ) /X i ns 
5710 Qlo=(QloUQlo2)/2 

5720 Qli1-(Kins»Apl^(Tinsi1-Tinsi2))/Xins 
5730 01i2=(Kins-»Apl ^(Tinsi2“Tinsi3))/Xins 
5740 0li=(Qli1+Qli2)/2 

5750 Rr = 2.0>( ( l-£cu)/(ApHEcu) )■♦“( 1 .0/Apl) 

5760 Qr = ( Sigma* (ABS( ( Tuo+273 . 1 5 ) * 4- ( Tu i *273 . 1 5 ) * 4 ) ) ) /Rr 

5770 Qdel -Qp-Qa I r-Qlol -Qr-Ql i 1 

5780! 

5790 ! CALCULATE DEAN NUMBER. 

5800! 

5810 De= .5*Red*( (Dc/2)/Ri ) * .5 

5820 PRINT USING “ 1 2X , " “ Viscos i ty Based on Tblk (Mu) - "“,2Z.3DE."“ (kg/m.s) 
" " " ;Mu 

5830 PRINT USING “12X, ""Therm Cond Based on Tblk (Kair) - "“,Z.3DE."“ (W/m.K)" 
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; K a I r 



S8^0 


PRINT 






5850 


PRINT USING '•l2X.""P3tm 




-••.6D.DD,-- (N/m*2)--’‘;Patm 


5860 


PRINT USING •'12X, •’••Pdel 




“\4D.DD/'- (N/m*2)-”-:Pdel 


5870 


PRINT USING ■•12X,""P1 




-■•,6D.DD/'- (N/m*2)“-“;Pl 


5880 


PRINT USING "12X, -“Mdot 




-••,Z.4D,“- (kg/s)--“;Mdot 


5890 


PRINT USING ‘M2X,"“Repipe 




-••,6D,D“:Repipe 


5900 


PRINT USING "12X, •‘••Red 




-••,6D.D-;Red 


5910 

5920! 


PRINT USING "12X/*^‘Rehd 




'• \6D.D*‘:Rehd 


5930 


! ONLY PRINT DEAN NUMBER 


WHEN HEATING CURVED SECTION. 


59^0! 








5950 


IF Itype»l THEN PRINT USING 


"12X.“‘De - “•.4D.D’:De 


5960 


PRINT 






5970 


PRINT USING ■•12X,'‘-Qp 


M N 


.3D.3D,'" (Watts)"",7X. ■••Oair - “•.3D.3D,"" (W 


a t ts ) 


•' "“ ; Op , Oa i r 






5980 


PRINT USING -12X,-’‘01ol • 




".2D. 3D,"" (Watts)*",7X.""Qlil ■ "".2D. 3D,” 


(Watts) -•••' :Qlo1 .Ql i1 






5990 


PRINT USING “12X,“'‘01o2 - 




".2D.3D,"" (Watts)"", 7X, ""01 i2 - "".2D. 3D,"" 


(Watts) :0lo2.01i2 






6000 


PRINT USING -12X,--Qlo - 


m 


".2D.3D,"" (Watts)"". 7X. ""01 1 - "".2D. 3D,"" 


( Wa t ts ) ; Q 1 o , Q 1 1 






6010 


PRINT USING “l2X,'”0r 


t* 


".2D.3D,"" (Watts)"", 7X,""0del - "".2D. 3D,"" 


( Wa t ts ) “ * " ; Or ,0del 






6020 


PRINT 






6030 


PRINT USING “12X.”Hin - 


m M 


.3D.DD,"" (W/m‘2C)"".7X.""Nuin - "".3D.DD,";H 


i n . Nui n 






6040 


PRINT USING '‘12X.--Hmid - 


m II 


.3D.DD,"" (W/m*2C)"".7X,""Numid - "",3D.DD,";H 


mid, f4umid 






6050 


PRINT USING “12X,""Hout - 


m m 


.3D.DD,"" (W/n>-2C)"",7X,""Nuout - "",3D.DD,"lH 


ou t , Nuou t 






6060 


PRINT USING “12X,""Havg - 


m m 


.3D.DD,"" (W/m'2C)"".7X. ""Nuavg - “".3D.DD,";H 


avg , Nuavg 






6070 


PRINT 






6080 


PRINT 






6090 


BEEP 






6100 


INPUT ••WILL THERE BE ANOTHER 


RUN? ( 1-Y.0'N)",Go_on 


6110 


IF Go onOO THEN 




6120 


PRINT USING 






6130 


IF lm=0 THEN WAIT 300 






61^0 


GOTO 1700 






6150 


END IF 






6160 ' 


PRINT USING “SX/’-END OF RUN 


HUM 


61 70 


PRINT USING 






6180 


PRINTER IS 1 






6190 


ASSIGN -iFiie TO ^ 






6200 


END 






6210! 








6220 


! CONCERTS EMF TO DEGREES 


CELSIUS. 


6230! 

62^0 


SUB Tvsv(U,T) 






6250 


COM /Co/ D( 7) ,Aa(76.2) 






6260 


Sum = 0 






6270 


FOR 1=0 TO 7 






6280 


Sum=Sum + D( I ) I 






6290 


NEXT I 






6300 


T = Sum 






6310 


SUBEND 






6320! 

6330 


! CALIBRATES T/C READINGS 
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6340! 

6350 DEF FNTcorr(T,I) 

6360 COM /Co/ D(7) ,Aa(76,2) 
6370 Tc=Aa(I.O) 

6380 FOR J=1 TO 2 
6390 Tc»Tc+Aa(I.J)*T‘J 
6400 NEXT J 
64 to RETURN Tc 
6420 FNEND 
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APPENDIX B 

SAMPLE CALCULATION 



Figure (B.l) shows the components of the heat balance in the test channel. A 
sample calculation follows that details the methods used in the data reduction program 
to estimate the value of the energy balance components and then calculate the 
hydraulic Reynolds number and average Nusselt number for each data set. A copy of 
an actual record of one set of experimental data for the concave test section precedes 
the sample calculations together with a list of the parameters necessary to make the 
sample calculation. 




Figure B.l Test Section Energy Balance. 
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1. SAMPLE CALCULATION DATA 



Month, Day, and Time: 09:16:15:14:5^ 

THE FOLLOWING DATA WERE RECORDED: 

Heater ^j^oltage (Vh) - 91.755 (V) 

Resistor Voltage (Vr) = 7.704 (V*) 

Precision Resistor (Rpr^ = 2.01729 lOhms) 

OUTER PLATE: 



T/C NufTiDer: 


: 1 


2 


'!{ 




5 


6 


T enp ( C ) : 


72.33 


69.0 7 


7 6 . LI 0 


68.01 


76.30 


72.90 


T/C Numoer: 


1 

/ 


6 


9 


10 


1 I 


12 


T emp ( C ) : 


79.99 


72. 13 


73.77 


74.26 


77.03 


75.00 


T/C Number: 


: 13 


14 


15 


16 


17 


1 8 


T emp ( C ) : 


74.32 


75.52 


75.70 


63.57 


62.45 


78.73 


T/C Number: 


13 


20 


21 




23 


24 


T emp ( C ) : 


79.13 


84.26 


85.87 


69^58 


76.27 


78.07 


T/C Number: 


25 


26 


27 


28 


29 


30 


T emp ( C ) : 


7?\S3 


80.33 


47.70 


79.25 


80.06 


83.47 


OUTER INSULATION: 












T/C Number: 


31 ( T 1 nso 1 ) 


32 ( T 1 nso2 ) 33 


( T 1 nso3 ) 




T emp ( C) : 


82, 


.38 


56.88 




40.01 




ORIFICE TEMP 


( To rf ) ” 


30.25 


(C) 








INLET TEMPERATURE: 












T/C Number: 


35 


36 


37 


38 






T emp ( C) : 


22.09 


22.07 


22.08 


22.07 






INNER PLATE: 














T/C Number: 


41 


42 


43 


44 


45 


46 


T emp ( C ) : 


23.89 


23.69 


23.40 


23.34 


23.21 


23.49 


T/C Number: 


47 


48 


49 


50 


51 


52 


T emp ( C) : 


23.28 


26.29 


23.56 


23.27 


24.03 


24.23 


T/C Number: 


53 


54 


55 


56 


57 


58 


T emp ( C) : 


23.71 


23.82 


23.59 


23.63 


23.56 


24.02 


T/C Number: 


59 


60 


61 


62 


63 


64 


T emp ( C) : 


23.80 


24.49 


24.27 


24.87 


24.68 


25.87 


T/C Number: 


65 


66 


67 


68 


69 


70 


T emp ( C) : 


25.56 


25.32 


25.16 


25.04 


24.95 


24.93 


INNER INSULATION: 












T/C Number: 


71 ( T 1 ns 1 1 ) 


72 (Tinsi2) 73 


( T 1 ns 1 3 ) 




Temp ( C) : 


24, 


,40 


24.13 




23.70 




OUTLET TEMPERATURE: 












T/C Number: 


74 


75 


76 


77 






T emp ( C ) : 


30.66 


30.73 


30.65 


30.69 
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THE FOLLOWING TEMPERATURES WERE CALCULATED: 



Average Outer Wall Temperature (Tuo) 

Average Inner Wail Temperature (Tui) 

Average Outlet Temperature (Tout) 

Average Inlet Temperature (Tin) 

Channel Inlet and Outlet Temp Diff (Tdiff) 
Average Bulk Temperature (TbU) 

Mean Temperature Difference (Tdel) 

Outer and Inner Wall Temp Difference (Tudiff) 
Average Local Channel Temp at Inlet (Tcin) 
Average Local Channel Temp at Middle (Tcmid) 
Average Local Channel Temp at Outlet (Tcout) - 
Avg Local Heated Wall Temp at Inlet (Tuin) 

Avg Local Heated Wall Temp Near Middle (Tumid) = 
Avg Local Heated Wall Temp at Outlet (Tuout) 



78.^6 ( C ) 
24. 16 (C) 
30.68 (C) 

22.08 (C) 
3.61 (C) 

26.38 (C) 

50. 08 (C) 

52.30 (C) 
22.44 (C) 

25.30 (C) 
30.32 (C) 
73.62 (C) 
73.06 (C) 
79.78 (C) 



THE FOLLOWING ORIFICE DATA WERE ENTERED: 



Pa tm ( 1 nHg ) 
29.86 



DPM( inH20) 
17.60 



PlM(inH20) REt 

16.90 1S000.0 



Do r f ( m ) A ( m * 2 ) Beta 

0.027318 S.361E-04 0.S277 



THE FOLLOWING DATA WERE CALCULATED: 



Orifice Expansion Factor (Y) 
Orifice Flou Coefficient (K) 
Density Based on Torf (Rho) 
Uiscosity Based on Torf (Mu) 
Viscosity Based on Tblk (Mu) 
Therm Cond Based on Tblk (Kair) 



0.98S9 

0.6326 

1.1099 (kg/m‘3) 

18.6I6E~06 (kg/m.s) 
13.434E-06 (kq/m.3) 
2.620E-02 (W/m.K) 



Pa tm 




- 1008A8.13 (N/m'2) 








Pde 1 




A379.58 (N/ra-2) 








PI 




= 966^2.74 (N/m*2) 








Mdo t 




= 0.8360 (ka/s) 








Rep ipe 


= 47611.8 








Red 




7697.2 








Rehd 




= 15013.9 








De 




397.9 








Qp 




350.411 (Watts) 


Qa 1 r = 


31 1 .991 


(Watts) 


Qlol 


= 


5.908 (Watts) 


Qlil = 


.063 


(Watts) 


01o2 


= 


3.908 (Watts) 


Qli2 = 


. 100 


(Watts) 


Olo 


= 


4.308 (Watts) 


fill = 


.001 


1 W .3 r t a i 


Or 


= 


1 . 396 (Watts) 


Odel = 


30 . 4SJ 


( i3a t ts ) 


H i n 




78.74 (W/m'2C) 


Nu 1 n 


= 37.67 




Hm 1 d 




84.3.9 (W/m'2C) 


Hum i d 


^ 4 0 . 0 3 




Hou t 


= 


81 .43 (W/m‘2C) 


f'luou t 


= 38.03 




Havg 


= 


80.47 ( W/m ■ 2C ) 


Nuavg 


= 38. OS 





END OF RUN 
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= 


91.755 V 


^pr 


= 


7.704 V 


V 


= 


2.01729 n 


^atm 


= 


29.86 in. Hg 


AP 


= 


17.60 in. H^O 


Pi 


= 


16.90 in. H^O 


T. 

in 


= 


22.08 °C 


T 

out 


= 


30.68 °C 


T 

wo 


= 


76.46 °C 


T . 

Wl 


= 


24.16 °C 


T. 

insil 


= 


24.40 °C 


T. . 

insi2 


= 


24.13 °C 


^insi3 


= 


23.70 °C 


T. 

insol 




82.38 °C 


T 

inso2 


= 


56.88 °C 


7. 

inso3 


= 


U 

o 

o 


Torf 




30.25 °C 


Kair 


= 


.0262 W m-°K 




= 


1007 J Kg-°K 


f'air 


= 


4.6971 X 10'8 > 


K. 

ins 


= 


.0380 W/m-°C 


ins 


= 


.0127 m 


^cu 


= 


0.12 


(T 


= 


5.669 X 10‘S W, 


P 


= 


.5277 


Y 


= 


1.40 






1.0 Kg-m, N-s^ 


R 


= 


286.98 Nm Kg- 


F 

WO-Wl 


= 


1.0 


■'^pl 


= 


.07742 nr 




=: 


.0016 nr 


A . 

pipe 


= 


.00211 nr 


D 

C 


= 


.00635 m 


D . 

pipe 


= 


.0518 m 


D 

o 




.027318 m 



wet 



= .5207 m 
= .297 m 



A- 5 



air 
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TEMPERATURE CALCULATIONS 



a. Average Temperature of the Outer (Inner) Wall 
^ Outer (Inner) Plate Temperatures 
wo wi Number of Outer (Inner) Plate Thermocouples 



= 76.46°C (24.16°C) 



b. Average Channel Outlet (Inlet) Temperature 
^ Outlet (Inlet) Temperatures 
out^ in^ Number of Outlet (Inlet) Thermocouples 



= 30.68°C (22.08°C) 



c. Channel Inlet to Outlet Temperature Difference 

T^ fr = T „ - T. = 30.68°C - 22.08°C = 8.6°C 

dill out in 

d. Bulk Temperature of the Flow Through the Channel 

T + T , 22.08X + 30.68X 

T. 2llt_ = = 26.38°C 

2.0 2.0 

e. Average Temperature Difference of the Heated Wall and Fluid Bulk Temperature 

" \o ■ \lk " ■'6-46°C - 26.38°C = o0.08°C 



f Temperature Difference Between the Outer and Inner Walls 
T . «• = T - T . = 76.46X - 24.16°C = 52.30°C 

wdill wo wi 



3. POWER CALCULATIONS 



Qn = 



a. Power supplied by the Heated Plate 

V, XV 91.755V X 7.704V 
= ^ = = 350.41 Watts 



R_ 



pr 



2.01729L2 
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b. Po>ver Loss Through the Outer Wall Insulation 
Qlc - <''2>(Ql,l + <3 .o2> 



- (l/2)[( 



T - T. 

insol inso2 \u \ 

AX. P> 



T. - T. 

1 / inso2 inso3 m. a i = 

^ AX ^ ins^^plJ 

ins 

^ (0.038 W, m"C) 

{l/2)(.07742m^) [(82.38 - 56.88(°C)) 

.0127m 



(56.88 - 40.01(°C))] = 4.908W 



c. Power Loss Through the Inner Wall Insulation 

Q,i = (l'2XQ,i, + Q,,) 

= (l/2)[(4j!!sil-lIumL)k A, 



AX., 



ins pi 



( 



A„,l - 



AX. P‘ 

ins 



, (0.038 W'm"C) 

(l/2)(.07742m-) - [(24.40 - 24.13(X)) + 24.13 - 23.70(°C))] = .081W 

.0127m 

d. Heat Loss by Radiation 



1. Radiation Resistance 



R = 2X( 



1 - 



Apl X Ccu 
I 2 

‘^pl ^CU 



•) + 



1 



A I X p 
pi wo W1 

1 



1) = -,( — -1) 

0.07742m- 0.12 



= 202.36 m 



2. Heat Radiated 
O - 'XT'vo'' - 

K 

(5.669 X 10 ^)[( 76.46 -t- 273)-^ - (24.16+ 273)-^] 



202.36 



1.996 W 



4. MASS FLOW RATE CALCULATIONS 



a. Pressure Conversions 

P , = 29.865 in Hg x 3376.8 = 100848.13 N 

AP = 17.60 in H 2 O x 248.84 = 4379.58 N nr 
P, = (16.90 in H,0 x 248.84) - 100848.13 = 96642.74 N 

b. Density of Air 

P, (96642.74) , 

p . s — ^ = 1.1099 Kg,m^ 

RT f (286.98) (30.25 + 273) 



c. E.xpansion Factor 

o4 

Y = 1 - (0.41 + 0.35p^) 

YPj 

, 4379.58 

= 1 - [0.41 + 0.35(0.5277)^] = 0.9859 

(1.401(96642.74) 



d. Mass Flow Rate 

m = YKA v/2g^P~^= (.9859)K(.0005861m^) 

X V2(l Kg-m,Xs-)( 1.1099 Kg nr)(4379.58 N'm‘) = .05669K Kg,s 
Iterate to find ni: 

Assume a hydraulic Reynolds number = 15000 



Re . 

pipe 



D. ,A . 

hd pipe 



) = ( 



.0518 X .0016 
.01229 X .00211 



P-'^pipe 



15000(3.1961) = 47942 



3.1961 X 



Interpolate in the appropriate table in Reference 15 for K, the (low coefiicient: 
K = .6329 



{Xote then the computer program used to reduce the data used an approximating 
polynomial to find the correct value of the flow coejflicient, K.) 



-r m = (.05669)(.6329) = .0359 Kg s 
Solve for the improved value of Re . ; 

^ pipe 

liiD. (.0359K.0518) 

^ T7- = 47315 

-\ipePair (-0021 1)( 18.616 X 10 <’) 

Iterate until successive values of m dilTer by less than .001: rh = .0360 
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5. REYNOLDS NUMBER CALCULATIONS 



a. Re Based on Pipe Diameter 

mD. (,0360)(.0518) 

Re . = Eiee. = H ’ = 47475 

A. u. (.00211X18.616X 10 ^) 

b. Re Based on Channel Height 



rh D 

Re, = ^ = 

A R . 

c^air 



(.0360)(.00635) 

(.0016X18.434 X io'6) 



= 7751 



c. Re Based on Hydraulic Diameter of the Channel 

(.0360X.01229) 



rh D, . 

Re = = 

^c^^air (.0016X18.434 X10'^) 



= 15000.8 



6. HEAT TRANSFERRED TO THE WORKING FLUID 

Qair = (Tout ■ Tin) = ('0360 Kg,s)(1007 J 'Kg-°C)(30.86 - 22.08)°C = 312 W 



7. 

h = 



AVERAGE HEAT TRANSFER COEFFICIENT CALCULATION 



Qair 

Apl^T 



312 W 



(.07742m")(50.08°C) 



= 80.5 W;m^-°C 



8. AVERAGE NUSSELT NUMBER CALCULATION 

hD.. (80.5 W/m^-“C)(.0 1229m-) 

^ ^ = 37.8 

K,.. .0262 W,m-’C 



9. DEAN NUMBER CALCULATION 
Re 



De = — ^ = 
2 



'Dc'2 

r. 



7751 

"1 



y 



.00635.2 

J97 



400.7 
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APPENDIX C 

EXPERIMENTAL UNCERTAINTY 



The experimental uncertainty in calculated parameters was deterntined in 
accordance with the method described by Kline and McClintock [Ref. 23]. 

The qualitative results shown here correspond to the data set and calculation of 
Appendix B. There was a significant decrease in the experimental uncertainty of the 
average Nusselt number for Reynolds numbers above 17,000 in this experimentation 
compared to that of Galyo [Ref 14: pp. 104-105] because of the changes made to the 
test channel that promoted mixing of the working fluid at the exit of the channel. 

The equations used to determine the experimental uncertainty of calculated 
parameters are shown below. 



I. DENSITY OF THE WORKING FLUID 



6p 6R 

— = [( — 
P R 



)* + ( 



6P 



1 ^2 



)" + ( 



6T 



orf .0.5 



r r 



1 



orf 



2. MASS FLOW RATE OF THE WORKING FLUID 

6m 6Y , 6K , 6A , 

— =[(( — )- + ( — )“ + ( —)- ) 
m \ k A 

6p . , 6AP n ^ 

+ (I + ( TIT* !l 

l>a„ 



HEAT TRANSFERRED TO THE WORKING FLUID 



6Q . 

^a.r 

Q • 



6 m 



6c_ 



= [( — )- + ( 

m Cp 

(Tout - Tin) 






)2 j 0.5 



6S 



4. AVERAGE HEAT TRANSFER COEFFICIENT 
6h 

T~ 



60 . , 5A . , 6AT , n , 

[( Ta»r )2 + / eL )2 + ( )2 J 0.5 

Q . ’ A . AT 

-<q - 



^air 



pi 



5. AVERAGE NUSSELT NUMBER 

6D 



6Ntx . 6h , 

= [( — )^ + ( 

>^\d h D,, 



bA ^2 ^ ^air ^2 j 0.5 



Kair 



6. HYDRAULIC REYNOLDS NUMBER 

6|i . 



6Re 6m , 6D. . , 

= [( — )- + ( 9 + ( 



6A 



Re 



m 



D 



hd 






+ (— S-) 



c \2 1 0.5 



Several parameters used in the calculations described herein had uncertainties 
that were constant over the entire range of testing. These quantities and their relative 
uncertainties follow. 



Quantity (B) 



Uncertainty ( ± 6B B) 



A 

A. 



pi 



P 

D, 



D 

K 



hd 



^^air 



atm 



R 

W. 



id 



Y 



air 



.0020 

.0408 

.0006 

.0020 

.0400 

.0391 

.0030 

.0038 

.0010 

.0003 

.0080 

.0020 

.0054 
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The approximating polynomials used to resolve thermocouple voltages into 
temperatures matched actual temperatures to within less than±0. 1'C over a range of 
temperatures from 20'’C to 85'’C; however, in some cases this temperature range was 
exceeded. A conservative estimate of the uncertainty of temperatures measured by 
thermocouples was set at ±0.rc, the accuracy of the thermocouples specified by the 
manufacturer of the thermocouple wire. The following uncertainties were used in 
calculations: 

Quantity (T) Uncertainty ( ± T) 



For the data taken at a hydraulic Reynolds number of 15000 in the concave test 
section the following additional uncertainties applied: 




0.1 ”C 
0.1 X 
0.1 ”C 
0.1 °C 
0.1 °C 
0.2 “C 
0.2 ”C 



AT 



Quantity (P) 



Uncertainty ( ± P) 



AP 



0.05 

0.10 
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The results of the uncertainty calculations at a Reynolds number of 15000 in the 
concave test section were; 



Quantity (B) 



Uncertainty ( ±6B/B) 



p 


.00445 


m 


.00544 


Q 

^air 


.02397 


h 


.02431 




.04620 




.05703 



The major contribution to experimental uncertainty in both the average Nusselt 
number and hydraulic Reynolds number was the uncertainty in cross sectional area of 
the channel. The hydraulic diameter of the channel was a function of channel cross 
sectional area, resulting in a relatively large uncertainty in the hydraulic diameter. 
These two parameters accounted for approximately fourty percent of the uncertainty in 
the average Nusselt number calculation, and about eighty percent of the uncertainty in 
the calculated hydraulic Reynolds number. 
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